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Abstract 

The relationship for Pulsed-Field-Gradient NMR between the amplitude I{t) 
of the spin echo and the molecular displacement X(t) is examined. I{t) of a sin- 
gle species in a simple solution is determined by the mean-square displacement 
X(t)^. With polydisperse species, or molecular probes in complex fluids show- 
ing memory effects, I{t) in general includes large contributions from all higher 
even moments Conditions under which the NMR signal is indeed de- 

termined by the molecular mean-square displacement are noted. A diagnostic 
that sometimes identifies when these conditions are not met is presented. 

1 Introduction 

Pulsed-Field-Gradient Nuclear Magnetic Resonance has proven to be an effective 
means of studying diffusion in simple and complex liquids. The quantity directly 
measured is the dependence of the amplitude of the spin echo on the strength 
and time separation of the gradient pulses. In simple non-viscous liquids, the 
strength of the spin echo corresponds, on a molecular level, to the mean-square 
molecular displacement and to the molecular self-diffusion coefficient. In com- 
plex fluids, matters change. The objective of this paper is to treat these changes 
on a fundamental level, obtaining relationships between spin echoes and the full 
probability distribution for single-molecule displacements. 

The starting point is the relationship between the intensity I{t) of the spin 
echo, the experimental parameters, and the molecular properties [I]. The ex- 
periment applies short-duration magnetic field gradients G separated by a time 
interval A, gradients enduring for a duration S. The significant time t is related 
to A and 6 by t = A — (5/3. The distribution function for molecular single- 
particle displacements X during time t is P(X, t). The relaxation of I(t) with 
increasing t is related to molecular motions hy^ 

I{t)/I{0)= y dXP(X,t)exp(27riq- X), (1) 

the wave vector being 

q=Gj5/i2n). (2) 
where molecules of interest have gyromagnetic ratio 7. 
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Callaghan and Coy[3] provide a slightly different but equivalent expression 
for the time dependence of I{t), which may be written 

g^''Hq,t) = ^^^exp[27r*q. X„(t)]^ , (3) 

in which g'^^'^'^ {q,t) is the self part of the dynamic structure factor, and in which 
X„ (t) is the displacement of the n*'' of the N particles during the time interval 
t. 

I{t) and g^^^\q,t) are representations of the incoherent (single-particle) dy- 
namic structure factor. In somewhat different systems with different time and 
distance scales, at least four experimental methodologies are in use to determine 
g^^^'iq^t). In addition to magnetic resonance methods, 5(1^) is obtained using 
(i) inelastic neutron scattering from dilute isotopically-tagged tracer polymers 
in a polymeric fluid [4j, (ii) optical tracer diffusion in which light scattering spec- 
troscopy is applied to dilute solutions of an intensely-scattering tracer polymer, 
e.g., polystyrene, in solutions of an isorefractive but possibly non-dilute matrix 
polymer, e.g. polyvinylmethylether:orthofluorotolueneL5, , and (iii) optical probe 
diffusion in which light scattering spectroscopy is used to observe intensely- 
scattering probe particles such as polystyrene spheres moving through solutions 
of weakly-scattering probe polymers^. While these methods have very different 
physical bases, they all measure the same underlying correlation function, seen 
on the right- hand-sides of eqns. [1] and [3l Theoretical considerations 7 applied 
to interpret relaxation in any of these methods therefore remain applicable for 
the other methods. (Note, however, that in treatments of neutron and light 
scattering the factors of 2tt in eqns. [2] and [3] are cancelled against each other; 
the wave vector of eqn. [5] would in light scattering be termed the inverse wave 
length.) 



2 Results 

We now consider relationships between g^^'^^Qjt) and P(X,i) for progressively 
more complicated systems. Conditions under which g^^^\q,t) is determined by 
the mean-square displacement are identified. A diagnostic that sometimes iden- 
tifies when g^^'^^iqjt) is not being determined by the mean-square displacement 
is presented. 

The simplest ideal system involves a solution of identical dilute particles 
whose motions are described by the Langevin equation. For this ideal case, the 
distribution of particle displacements has a Gaussian form 

P(X,0= (y((X)2)) '^'exp(-3(X)V2((X)2)), (4) 

with the brackets (• • ■ ) indicating an ensemble average. The mean-square par- 
ticle displacements satisfy ((X)'^) = 6Dt. 
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For identical Langevin particles, the displacements are also a Markoff pro- 
cess, meaning that the particle displacements during non-overlapping time in- 
tervals are independent from each other. The Gaussian and MarkofF properties 
are independent of each other: A statistical process may be Gaussian but not 
Markofhan or may be Markoffian but not Gaussian. For the ideal system, Berne 
and Pecora[B] evaluate eq. [T]and compute g^^^^q,t), finding 



The relaxation of the incoherent structure factor is determined by the mean- 
square particle displacements. 

The extremely restrictive nature of the conditions under which eq. [S] is valid 
is often not emphasized. Figures [T] and [5] demonstrate the complete failure 
of eq. [5] when the conditions have been very slightly relaxed. Figure [1] shows 
f^^^"-* (q, t) computed for a nominal bidisperse system of Brownian particles. Fig- 
ureHplots (solid line) - ln(g(i'') (g, <))/(47r2g2 

), which would be the mean-square 
displacement if eq. [5]were correct. Figure [5] compares this form with the true 
mean-square displacement of the bidisperse particles (dashed line). If eq.[5]were 
correct for this system, the two lines in the Figure would be exactly the same. 
In fact, except that they share a common origin and initial slope, the lines are 
entirely different, because the conditions for the validity of eq.[S]are not satisfied 
by bidisperse Langevin particles. Particles of each species retain their sizes at 
all times, so AX(t) is not a Markoff process. Larger and smaller particles each 
persist in diffusing more rapidly or less rapidly, no matter for how long they 
are observed. In terms of eq. HI for each particle P(X, t) is a Gaussian, but the 
average of -P(X, t) over all particles is not a Gaussian, but is instead the sum of 
two Gaussians. 

The relationship between particle displacements, Gaussian processes, Markoff 
processes, and invocations of the Central Limit Theorem is made precise by 
Doob's theorem[9] and the discussion around it. This matter has previously 
been analyzed in depth for light scattering spectroscopy [7], so key issues are 
only summarized here. Note, however, that light scattering spectroscopy in- 
volves coherent scattering; correlations involving two different molecules at two 
different times contribute to the light scattering spectrum. Magnetic resonance 
involves the incoherent structure factor; terms in the light scattering treatment 
of (?^^^ (q, t) that involve two different molecules average in the magnetic reso- 
nance treatment exactly to zero. With this reminder, only a few key lines of the 
analysis of light scattering spectroscopy need to be repeated here. 

Doob originally considered a random process u{t), a random process being 
a time-dependent variable whose steps from time to time were random but had 
statistical properties. u(t) at a given time had a Gaussian random distribution, 
and its distribution at two times had a joint Gaussian distribution. An example 
of u{t) is the q*^ spatial Fourier component of the single-particle density 



(5) 




(6) 



n 
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Here R„(t) is the location of particle n at time t. The remaining discussion is 
phrased in terms of aq rather than u. 

Because the molecular locations are independently randomly distributed, all 
with the same distribution, and we are only considering the incoherent structure 
factor, the Central Limit Theorem guarantees that aq has a Gaussian random 
distribution. For the incoherent structure factor, the Central Limit Theorem 
also shows that the displacements X„(t) = R„(t) — Rn(0) are independent, all 
having the same distribution, so the incoherent part of aq (t) at any pair of times 
has a joint random distribution. 

In the case that aq{t) is also a MarkofF process, meaning that the probability 
distribution function for aq{t) at each of a series of times depends only on the 
value of aq{t) at the immediately prior time, Doob demonstrates that 

(a,(0)a,(t)) = (|a,(0)p)exp(-rK|), (7) 

except for the special case {aq{0)aq{t)) — 0. The demonstration applies equally 
to the coherent and incoherent parts of {aq{Q)aq{t)). For a Gaussian MarkofF 
process, the relaxation is a simple exponential in t. 

The Central Limit Theorem does not guarantee that a process has Markoff 
behavior, because the displacement of each particle during each time interval 
may be correlated with the displacement of the same particle during the prior 
time interval. The absence of MarkofF behavior is seen in velocimetry studies, 
in which the displacements of tracer particles during sequential time intervals 
are very heavily correlated and, correspondingly, g^^^\q,t) shows not an expo- 
nential relaxation but instead, as it happens in this case, sinusoidal oscillations. 
Contrapositively, if the relaxation of 5^^"' {q, t) is not a simple exponential, then 
the fluctuations in aq{t) are not described by a Gaussian Markoff process. The 
lack of a simple exponential relaxation for g'^^'^^q^t) is an infallible diagnostic 
that the assumptions leading to eq. [5] are not applicable in the system under 
consideration. 

In systems in which g'^^'^\q,t) — exp{4'K^ q^ Dt) , there is no need to take 
short-time or small-g limits, because the slopes of log((7*^^*' (g, t)) as functions of 
t or q"^ are the same at all t and q. Contrariwise, if experimentally \og{g^^'^^ {q, t)) 
is not linear in t, Doob's theorem guarantees that the derivation leading to eq. 
[5] is not applicable to the system under study. 

What happens in the non-ideal cases in which Doob's theorem is not ap- 
plicable? Magnetic resonance and scattering both can measure g^^^^ {q, t) as a 
function of g or as a function of t. In real applications, practical differences 
arise. Light scattering experiments observe light scattered through a single an- 
gle; g^^'^\q,t) is automatically simultaneously computed at a large number of 
different times. Observing g^^^'^q, t) at several q with light scattering requires a 
disproportionate increase in experimental effort. A single pulsed field gradient 
NMR cycle observes g^^'^^{q, t) for a single pair of values g, t. NMR studies can 
vary t, and can vary q by changing G. Leaving t fixed avoids complications 
attendant to various relaxations |:2_, encouraging NMR studies to vary q rather 
than t. 
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For small values of its exponential's arguments, and noting that symmetry 
eliminates terms odd in X, a Taylor series expansion of the exponential in eq. 
[3] gives 



4! 



7V-i^((q.X,(<))4)+0(^ 



the lead terms of which are 



(8) 



(9) 



Here 



N 



((q-x,wr 



(10) 



Equation IH] becomes accurate asymptotically in the limit of short time or in 
the limit of small wavevector. The short time limit is approached by evaluating 
limt^o dg'^^\q, t)/dt, as by fitting g^^^{q, t) to a cumulant expansion [TUl [TT| . The 
small-g limit is approached by going to small angle or small field gradient. 

If actually attained, the small-t and small-g limits are physically very dif- 
ferent. In the true small-time limit, particles have not had an opportunity to 
move very far; all relaxation process are still ongoing. In the true small-q limit, 
particles must displace through large distances to contribute to the relaxation, 
which takes a long time, so relaxation processes and memory-function effects 
will all have decayed. 

It might be tempting to replace eqn. [9]with a supposed exponential approx- 
imant 

(2^)2 



/i^)(g,i)=exp 



-q\^X{t)f 



(11) 



However, a full calculation shows that eqn. [TT]is incorrect in its first non-trivial 
term, namely its Taylor series expansion is already wrong in its q^ term. The er- 



ror is revealed by noting that all quantities in exp(- 



-q X{tY) are constants 



that may be factored out of the ensemble average of eqn. [3l yielding 

(2^)2 



exp(- 



N 



7V-i^ <^exp(27r2q- Xj (<) 



-q'X{tf) 



q'xm 



(12) 



A Taylor series expansion of the second exponential leads to 
g^'^\q,t)^ [exp(- 
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(13) 



The first correction appears at the level. The exponential form is no more 
accurate that the linear form it was proposed to replace. 
Iteration of the exponentiation process leads to a series 
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-15X4X2+ 30X2 ) + ... . (14) 



Equation [T3] gives the relationship between the incoherent structure factor and 
various moments of the particle displacements. The equation reduces to eqn. 
[TT] if and only if P(X, t) is a Gaussian, because in that case X^ — 3X2 g^j-^^ 
corresponding terms of higher order in X all vanish. If P(X, t) is a sum of 
Gaussians, as is the case for a polydisperse system, or is not Gaussian at all, 
eqn [Ml does not reduce to eqn. [TTl 

Methods used to interpret light scattering spectra suggest two approaches 
to use equation 1141 to interpret intensity relaxations I{t). 

First, light scattering spectra, which are obtained at fixed q and a large 
number of delay times, are routinely interpreted by means of time cumulants 

K,. ^ lim(-,^)-"^::M|^ (15) 

Ki is the light-scattering-intensity- weighted average diffusion coefficient D, while 
K2 gives the mean-square range of diffusion coefficients — D . Similar iden- 
tifications can be made for the higher time cumulants. In practice, the Kn are 
determined by fitting g^^''^ (q, t) to a power series in t. The number of cumulants 
that can be determined accurately is determined by the range of times studied 
and the signal-to-noise ratio, but is never large. 

Second, by analogy with time cumulants, one may propose the use of wave- 
vector cumulants 

Q„ = (2n)!(^-j hm . (16) 

By direct calculation. 



Qi = X(t)2 (17) 
is the mean-square displacement at time while 

Q2 = (X4-3X2') (18) 

gives the deviation from a Gaussian form of the displacement distribution. If 
P(X, t) has a simple Gaussian form, Q2 and all higher wavve-vector cumulants 
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Qn, n > 2, vanish. By analogy with the time cumulants, the second wavevector 
cumulant is usefully normalized as 

(19) 

The third wavevector cumulant is 

Q3 = (X6 - 15X4 X2 + 30X2^). (20) 

Just as time cumulants are usefully obtained by fitting g*-^*-* (g, i) as a function 
of time to a power series in t, so also wavevector cumulants should usefully be 
accessible by fitting g^^\q,t) as a function of q to a power series in q^. Linear- 
least-square fits to power series will reliably determine the leading initial slopes 
of g^^^\q,t) as a function of t or q^, even when the higher-order cumulants 
(and, hence, the degree of curvature of \n{g'^^'^\q,t))) are considerably different 
from zero. Procedures for ascertaining the number of time cumulants that may 
actually be extracted from a given spectrumjl2 should remain applicable as a 
path to determining how many wavevector cumulants can be extracted from a 
given g'-^^\q,t). 

3 Discussion 

In the above, the relationship between the pulsed-field-gradient spin echo inten- 
sity and the displacement distribution function of the observed molecules was 
discussed. Only for monodisperse Brownian particles, for which (7^^^^(q,t) is 
a single exponential in q^ and t, is g^^^'^lq^t) determined by the mean-square 
particle displacement. If g^^*-* (q, t) is not a single exponential in t, g^'^'^^q^t) 
receives contributions from the higher moments X^", n > 2 of the displacement 
distribution. In general g''^^\q^t) depends on all higher even moments X^" of 
the distribution function for particle displacements. An expansion for (?^^*-'((7, t) 
in terms of the moments is presented. Analysis of g^^''-' ((7, t) in terms of time or 
wavevector cumulants was considered. 
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Figure 1 : Computed bimodal relaxation spectrum for a nominal mixture of two 
dilute Brownian species whose diffusion coefficients differ 30-fold. 
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Figure 2: Comparison of mean-square molecular displacement X[t)-^ with 
PFGNMR spin-echo amplitude plotted as — ln(7(t)) for the same mixture. 
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